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ABSTRACT 

In this paper we present the first results from implementing two scalar-tensor mod¬ 
ified gravity theories, the symmetron and the Hu-Sawicki /(I?)-gravity model, into 
a hydrodynamic N-body code with dark matter particles and a baryonic ideal gas. 
The study is a continuation of previous work where the symmetron and f(R) have 
been successfully implemented in the RAMSES code, but for dark matter only. By 
running simulations, we show that the deviation from ACDM in these models for the 
gas density profiles are significantly lower than the dark matter equivalents. When 
it comes to the matter power-spectrum we find that hydrodynamic simulations agree 
very well with dark matter only simulations as long as we consider scales larger than 
k ~ 0.5 h/Mpc. In general the effects of modified gravity on the baryonic gas is found 
to not always mirror the effects it has on the dark matter. The largest signature is 
found when considering temperature profiles. We find that the gas temperatures in 
the modified gravity model studied here show deviations, when compared to ACDM, 
that can be a factor of a few larger than the deviations found in density profiles and 
power spectra. 


1 INTRODUCTION 

The current standard model of cosmology is the ACDM 
model, and even though it is remarkably accurate in many 
aspects, it still has several problems that remains to be 
solved (Kroupa et al. 2010). Furthermore, ACDM has no 


proper motivation for adding the A constant other than the 
fact that it gives rise to the acceleration of the expansion 


of the Universe (Peebles & Ratra 20031. This can also be 


explained by introducing a new fluid known as dark en¬ 
ergy (Copeland, Sami & Tsujika wa|2006 1 to the matter con¬ 
tent of the Universe. Alternatively, the acceleration of the 
Universe might be a signal that gravity is modified on the 
largest scales. Modifying gravity is usually done by modi¬ 
fying the Einstein-Hilbert Lagrangian density £eh = R by 
replacing it with a more general function including terms 
of higher-order in derivatives of the metric (R 2 , R^R^ 1 ', 
Rap l i''R a ^ IJ ' v • ■ •) ° r by introducing new dynamical degrees 
of freedom, such as scalar fields (ICapozziello & de L aurentis| 
20111, coupled to the matter sector (Brookfi eld et al.|2006 |. 

General Relativity have been thoroughly tested in the 
lab (|Hoyle et al.||2004| |Dimopoulos et aT7||200T[) and in the 


solar-system and no deviation have so far been found (Will 


2014 Bertotti, less & Tortora 2003). This places strong 


constrains on any model that seeks to modify gravity and 
for such models to have any cosmological signatures, apart 
from modifying the background expansion history, a screen¬ 
ing mechanism is required. A screening mechanism ( |Khoury] 
2010 Brax et al.|2012 1 is a way of suppressing the effects of 


modifications of gravity in high density regions, compared to 
the critical density, such as on Earth and in the solar-system. 


In this paper, we will take a look at the symmetron (Hin- 
terbichler fe Khoury||2010 1 and /(A)-gravity models (Hu & 


Sawicki 2007), which are but two of numerous scalar-tensor 


held theories that possess such a screening mechanism. 

With the ever increasing number of theoretical mod¬ 
ified gravity models it is important to find ways to com¬ 
pare these models to observations in order to exclude the 
ones that are not viable. Local gravity experiments typically 
gives constrains on modified gravity models that translates 
into signatures being in the non-linear regime of structure 
formation. Therefore one of the ways to find useful observ¬ 
ables is to use N-body simulations. So far the majority of 
the cosmological N-body community that works with mod¬ 
ified gravity has focused solely on simulations made with 


Zhao| 2008| Li, Mota & Barrow|20iy Li et al. 

2012; Zhao, 

Li & Koyama 2011; Barreira et al. 2013 Puc 

lwein, Baldi 


& Springel 

2013; Schmidt||2009 

Brax et al.||2013| Li, Zhao 

& Koyama 

2013; Llinares & A 

ota 2014| Gronke, Llinares 


in doing so we neglect a wealth of physics, after all what 
is observed are photons emitted from luminous matter. To 
go beyond dark matter and try to describe the plethora of 
baryonic effects that take place in out Universe is very chal¬ 
lenging. Over the last decade, several codes have been writ¬ 
ten that combine N-body simulations with hydrodynamics 
( Teyssier|2002 Springel|2005 Cen 1992 Vogelsberger et al. 
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2 Amir Hammami et al. 


2014 Ryu et al. 19931 and include recipes for handling star 
formation and feedback processes. 

So far there has only been, to our knowledge, two at¬ 
tempts at combining modified gravity with hydrodynamics 
in N-body codes (jPuchwein, Baldi & Springel 2013; Arnold, 


Puchwein & Springel |2014 1. In this paper we show the first 

results from a hydrodynamic N-body code with the sym- 
metron model, and also reproduce and expand upon the 


same /(R)-gravity theory as presented in Puchwein, Baldi 
|fc Springel] ( |2013[ > and |Arnold, Puchwein &; Springel| ( |2014| ). 
The code that this paper is based on is a slight modification 


of ISIS (Llinares, Mota & Winther 2014), which in turn is 


a modification of the R AMSES hydrodynamic N-body code 


(Teyssier 2002). This paper will only concern itself of the 


modifications made in order to combine the modified grav¬ 
ity part of the N-body code with the hydrodynamic part 
of RAMSES. For more on the implementation of the scalar 
fields and other technicalities we refer the reader to |Llinares,| 
|Mota fe Winther] ( |2014[ ). The aim of this paper is to study 
the effects of modified gravity on the baryonic gas, compared 
to the effects that modified gravity has on the dark matter 
physics. For this purpose we will focus on simple observables 
such as the matter power spectrum and density and temper¬ 
ature profiles for both the dark matter and gas components 
separately. 

The structure of this paper is as follows. In Section [2] 
we give an introduction to scalar-tensor theories of gravity 
and take a look at how the equations differ from those of 
standard gravity. Section [3] briefly details the code imple¬ 
mentations and the run parameters of our simulations. We 
show our results for density profiles in Section[4j for temper¬ 
ature profiles in Section [5] and for power spectra in Section 
[6] We finish the paper with a short discussion in Section [7] 


2 SCALAR-TENSOR THEORIES OF GRAVITY 

We are interested in scalar-tensor theories that can are de¬ 
fined by the following action (Sotiriou 2006; Fujii & Maeda 
20031): 


S = J d 4 xW~9 


~^ M pi - ^d'tpdiip - V(tp) 


(1) 


T (c/[XL/ , \Ih ) ; 


where g is the determinant of the metric tensor, g in the 
Einstein frame, which is related to the Jordan frame metric 
tensor g by a conformal factor A(ip), 


9n» = A 2 (ip)g tJj 


( 2 ) 


and R is the Ricci scalar. The conformal factor satisfy A ~ 1 
for all the models we consider in this paper and we will use 
this approximation throughout. For more on these frames 
and the transformations between them and possible errors 
see|Faraoni, Gunzig fe Nardone| (|1999|) and |Brown fe Ham-| 


( 2012 ). 


Varying the action with respect to the scalar field and 
the metric gives us the scalar field equation of motion and 
the stress energy-tensor respectively, 


Dip = V'(ip) - A'(ip)T (m) , 


(3) 


where T ( "b i s the trace of the stress energy tensor, 


rp(m) 


= g^T^\ and 

T„u = AWT™ + 

= A(ip) [(P + p)ui_ l u„ + Pg^v] 

+ v^v.V' - Qstyaw + V(ip) 


(4) 


Note that the total stress-energy tensor is covariantly 
conserved, while the scalar field component itself is not 

VT^ ^ 0. 


We will now briefly go through the derivation of the 
fluid equations we get by using basic principles on the above 
action, and show what form they take for the symmetron 
model and the f(R)~ theory. 


2.1 The fluid equations 


We can derive the fluid equations of scalar-tensor theories 
from basic principles. We use the general action Eq. |lj and 
stress-energy tensor Eq. and start by computing the co¬ 
variant derivative of the stress-energy tensor, which is a con¬ 
served quantity (Misner, Thorne fe Wheeler|1973) 


= 0 . 


(5) 


We work in the Newtonian Gauge, 


ds 2 = —(1 + 2 &)dt 2 + a 2 (l — 2Q)8ijdx l dx ^, 


( 6 ) 


and from the conservation of the stress-energy tensor and 


by imposing the quasi-static limit (Llinares & Mota 2013 


Noller, von Braun-Bates & Ferreira 20141 for the scalar field, 


in which time-derivatives of the scalar field is ignored relative 
to spatial gradients, we find the fluid equations for scalar- 
tensor theories of gravity, 


+ V(vp) + 3Hp = 0, 
dv 


a 2 (P + P) [ Hv + ^ + 0 ' V)u + 

+VP + = o, 


d l +2 HE+- 
at p 


A(^) 

P 


(7) 

( 8 ) 


Vu = -(u-V)$-^M (u . V )V>. 


(9) 


We now perform a change of variables by implement¬ 
ing a variation of the so-called super-comoving coordinates, 


introduced by Martel & Shapiro (1998) 


ip = aip, 


P = a P, 


p = a p, 

<1 = a 2 ®, 


E = aE , 


( 10 ) 

( 11 ) 


where the tildes represent quantities in the super-comoving 
coordinates. The purpose of these coordinates is to eliminate 
unwanted dependencies on the scale and Hubble factors. The 
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equations then take the forn£] 

§+V(up) = 0, (12) 

§ + («-V)u = -ivP-V$-^MvVi, (13) 

dt p A(ip) 

^r+v- VE+^r- Vw =-(fl-V)4> - ^W-v- vjj, (14) 

dt p A(ip) 

which are the equations we have implemented in the N-body 
code. 


2.2 Symmetron model 


The symmetron model was introduced by (Hintcrbichler fe| 


Khoury (2010) as a new screening mechanism, similar to the 
one found in the Chameleon models (|Khoury fe Weltman| 


2004||Mota &: Shaw|2007). It utilises a screening mechanism 


so that we recover General Relativity in regions of high den¬ 
sity (such as the solar system where gravity has been very 
well tested) whereas we get an order one modification of 
gravity in low density regions. This is done by introducing 
a potential on the symmetry breaking form 


vw = Vo- ipV 2 + \\ip\ 


(15) 


where ip is the scalar held, p is a mass scale and A a dimen¬ 
sionless parameter. 

The coupling factor just mentioned is also chosen to be 
symmetric in the same manner as the potential 


AW = 1 + \ 



2 


where M is another mass scale. 

From the held equation we have that the dynamics of 
the held is determined by an effective potential given by, 

VeS {ip) = V 0 + 1 -[^-p 2 )ip 2 + 1 -\ip\ (16) 

In regions of high density ( p m M 2 p 2 ) the held is driven 
towards the minimum ip = 0, while in regions of low density 
we get a minimum at ipo = ±p^J^ for which the held will 
reside close to. The fifth-force (see below) is proportional to 
the local value of the scalar held so in high density regions 
ip fa 0 and it will be suppressed. 

We want, however, to work with slightly other param¬ 
eters, for which the physical interpretation is more clear, as 
presented in 

ing our free parameters p, M and A to /?, Ao and ossb like 


n Mpl'lj) 0 

; M 2 ’ 

(17) 

3 37/oH m M“; 

«SSB - MV , 

(18) 

A 2 - 1 

Ao “ V' 

(19) 


Winther, Mota & Li (2012). This entails chang- 


Now P represents the strength of the scalar hfth-force (rela¬ 
tive to the gravitational force), ossb is the expansion factor 


1 This results is reached by excluding terms of second order and 
assuming static pressure and fields. 


at the time of symmetry breaking, and is also related to 
the density at which the screening mechanism kicks in via 
the relation pssb = flmoPcoagSB! an d i s the ran g e of the 
scalar fifth-force in units of Mpc/h. Further the scalar field 
itself is replaced by a dimensionless scalar field \ by 

ip = ipox • (20) 


The equation of motion for this scalar field in the quasi-static 
limit is (Llinares, Mota & Winther]2014 1 


vr2 _ a 

X “ 2A 0 


r-rr 


Pm 3 

— + x - X 

Pm 


( 21 ) 


Simulations beyond the static limit were presented in 
|Llinares fc Mota| ( |2013||2014[ ), finding only sub-percent dif¬ 
ferences between the static and non-static solutions. 

The fifth-force in the symmetron model takes the form 


= 6 n^^^xVx, 

a SSB 


( 22 ) 


which is how it is represented in our code. For more on 
the symmetron model we refer to |Hinterbichler Khouryl 
( 2010 ). 


2.3 f(R)- gravity 


All f(R)- gravity theories revolve around promoting the 
Ricci scalar in the Einstein-Hilbert action to a function of 
the Ricci scalar instead. In the N-body code the Hu-Sawicki 
f(R) model (Hu & Sawic.ki 2007) has been implemented and 
is this is the model we focus on in this paper. The action for 
f{R)~ gravity is given by 


S = 




R + f(R) 

167 tG 


+ £n 


j4 

a x, 


and in the Hu-Sawicki model f(R) has the form 


f(R) = ~r 


, ci (R/m 2 ) n 
1 + C2(R/m. 2 ) r 


(23) 


(24) 


where n, ci and C 2 are the free parameters and m 2 = 
We can reduce the number of free parameters from 
three to two by demanding that Ci = 6 C 2 (to yield dark 


energy) as demonstrated in Hu & Sawicki (20071. Further it 
is convenient to use a parameter fuo instead of C 2 (LI inares,| 
Mota & Winther 2014) 

71 + 1 


/r 0 = — 


6nflA 

C2flm 


Ha 


(25) 


^ 3(H m -(- 4 Ha) . 

By applying the conformal transformation Eq. § to the 
action Eq. (231 using 


AW = e Mpl 


(26) 


with p = 1/V6 we recover the general scalar-tensor theory 
action Eq. 0- This model is then further characterised by 
the fu function 

f R = A 2 W-^*~- (27) 

Mpi 

In order to get the field-equation on a convenient form for 
a numerical implementation we must perform a change of 
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variables. As justified in Oyaizu (2008) we introduce a sub¬ 
stitution of the form 


r 2 u 

JR = ~a e , 


(28) 


combining this with Eq. (271 we see that this is essentially 
a substitution like 

l 2 A/pi 'll ,— , 2 A/pi Ur—7 

i/> = a ^-e => = a — e Vu. (29) 


In Oyaizu (2008) it has been shown that the equation of 
motion for this scalar held is then 


V • (e“V«) = f! m aHS 


-a* 1 + 4 


(/ 5-l)+(l + 4a 3 ^ 

S/a 


12 ?) 


a 


12?) 

2 j- \5ITT 


(30) 


'/so) ; 


where /ro is related to the range of the scalar Held via Ao oc 
1 /vTro ■ The fifth-force is given by 


= “/vv; 

^ A(iP) A/pi 

= ^ V “’ 


(31) 


which is how it is represented in our code. For more on f(R)~ 
theories of gravity see the review by|de Felice & Tsujikawa| 
( 2010 ). 


3 CODE IMPLEMENTATION AND RUN 
PARAMETERS 


Implementing scalar-tensor theories of gravity to the hydro- 
dynamic part of the N-body code is rather straightforward, 
thanks to the scalar-tensor theories all giving contributions 
as a fifth-force and the fact that RAMSES, which ISIS is 
based on, has been widely used, thoroughly tested and opti¬ 
mised. Wherever the code normally works with the gravita¬ 
tional force, call it Fgr, we simply replace it with an effective 
force F e g that includes the effects of modified gravity 

F e a = Fgr + Fiji, (32) 


which is then naturally permeated throughout the code. For 
more on the implementation of the scalar field solver itself 
see Llinares, Mota & Winther (2014). In the next two sub¬ 


sections we present our initial conditions and other param¬ 
eters we used in our runs. For both models we also made a 
run with standard gravity (ACDM) to use as a base of com¬ 
parison. All initial conditions have been generated by using 


the Grafic code, which is a part of COSMICS (Bertschinger 


19991, based on the parameters described below. 


3.1 Parameters for the symmetron simulations 

The symmetron simulations were run using 1024 cores, 256 3 
dark matter particles, with a box width of 256 Mpc/h and 
six levels of refinements. The background cosmology is a 
standard ACDM background with h = 0.65, I2a = 0.65, 
flm = 0.35 and fli, = 0.05. 

The symmetron parameters are presented in Table [I] 
These parameters were chosen to focus on various aspects 
of the symmetron models. Varying f) changes the strength of 


Table 1. Overview of the model parameters for the symmetron 
and f(R) models. 


Name 

P 

“SSB 


Sym_A 

1.0 

0.5 

1.0 

Sym_B 

1.0 

0.33 

1.0 

Sym_C 

2.0 

0.5 

1.0 

Sym_D 

1.0 

0.25 

1.0 

Name 

Pro 

n 


FofR04 

10 -4 

i 


FofR05 

10“ 5 

i 


FofR06 

10 -6 

i 



the fifth-force while changing assB changes the time of sym¬ 
metry breaking (i.e. it is the scale-factor for which the fifth- 
force kicks in) and also the density criteria for the screening 
mechanism. 


3.2 Parameters for the f(R) simulations 

The f(R) simulations were run using 1024 cores, 256 3 dark 
matter particles with a box width of 200 Mpc/h and eight 
levels of refinements. The background cosmology is a stan¬ 
dard ACDM background with h = 0.70, I2a = 0.727, 
12 m = 0.272 and I2 6 = 0.0450 

The f(R) parameters are presented in Table [l] These 
parameters as briefly introduced above were chosen as they 
give the full range of effects found in the model: from almost 
no screening and large deviations from ACDM for FofR04 
to much screening and small deviations from ACDM for 
FofR06. The main effect of changing f ro is to change the 
range over which the fifth-force is acting on and also the 
density threshold for screening. 


4 DENSITY PROFILES 


In this section, we present density profiles for multiple ha¬ 
los identified by using the the Rockstar code developed by 


Behroozi, Wechsler & Wu (2013) and incorporated into the 
YT-Project |Turk et al.| ( |2011[ ). We filter all the halos that 
have not reached a relaxed state. This is done by follow¬ 


ing the methods described in Neto et al. (20071 and Shaw 
|et al.| ( |2006| ), where we use relations between the kinetic- 
and potential energy and the surface pressure to determine 
if a halo is relaxed or not. We used the method presented 
by |Gronke( Llinares fe Mota| |20 14|) to take into account the 
effects of modified gravity in the virialization state of the 
halos. A halo is defined to be relaxed if 1 1T Jj Es < 0.2. See 
|Gro nke, Llinares & Mota] (j2014j for more on this limit and 
its implications. 

We are also interested in seeing how the behaviour de¬ 
pends on the mass of the halos. We will therefore study 
halos with masses in the ranges 1-5x10 13 /i _1 A/q and 1- 
5x 10 14 /t -1 A/q. The density profiles are calculated by bin¬ 
ning dark matter particles and the baryonic gas density in 


2 These ACDM parameters were chosen to coincide with those of 


Puchwein, Baldi & Springel 12013 \ to enable cross-checking. 
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Dark Matter Gas 



Figure 1. The top figures show the deviations from ACDM for our symmetron models for CDM and baryons respectively for all density 
profiles with a mass between 10 14 h~ 1 Mq and 5 X 10 14 /i _ 1 Mg. The bottom figures show the same for for all density profiles with a 
mass between 10 13 h~ x MQ and 5 X 10 13 /i _ 1 Mg. 


annular bins for each halo, then averaging over all halos. We 
focus only on the present day epoch which corresponds to 


z = 0. 


Our calculated density profiles are averages of all den¬ 
sity profiles of the proper size, ranging from 10% of the viri- 
alization radius, r = O.lffeooc, to ten times the virialization 
radius, r = 10f?200c- This range was chosen to properly catch 
all behaviours of the fifth-force on the dark matter and gas 
halos while also avoiding the inner regions of the halos where 
the resolution of our simulations is not sufficient. 


In Fig. [l] we present the deviations from ACDM of the 
density profiles for the dark matter and the baryonic gas for 
the symmetron model, while in Fig. [2] we present the same 
figures for the f(R )~gravity model. 

The dark matter density profiles for modified gravity 
show in general stronger clustering at the outskirts of the 
halos than the inner regions. The additional force from the 
scalar field will increase the rate at which dark matter and 
the gas collapses towards the centre. When the density cri¬ 
teria is then met and the fifth-force is screened the dark 
matter particles will have gained such velocities that they 
overshoot the centre of the halo and does not cluster there. 
The particles that overshoot the centre will then end up in 
the outskirts as we see in the density profile deviations. 


We observe that overall that the gas profiles behave 
closer to the ACDM counterpart then we find in the dark 
matter profiles. In previous works similar effects on the 


power spectra have been attributed to AGN feedback (Puch 
wein, Baldi & Springel 2013). In our simulations however we 


do not have these effects, as we chose to implement only the 
most rudimentary of baryonic components, a perfect fluid, 


clearly this effect is much more intrinsic to baryons than 
before assumed. 

What we see here might be an environmental effect, 
from the dark matter, on the gas component. The dark mat¬ 
ter collapses faster than the gas component, due to the CDM 
being collisionless while the gas components are hindered 
from collapsing due to friction and pressure as enforced by 
the Euler equations Eq. \12\ -Eq. ( | 13| ) , this means that the 
dark matter will cluster faster than the gas and reach higher 
densities earlier than the gas. The total density, dark mat¬ 
ter and gas, will trigger the density criteria that turns on 
the screening mechanism before the baryonic gas has had 
a chance to collapse as much as the dark matter. In other 
words, dark matter cuts off the fifth-force before it has had 
the time to work on the gas component to the same effect 
as it does on the dark matter. 

As seen in Table[l] we remember that all our symmetron 
simulations have the same coupling strength and force range, 
but with varying symmetry breaking criteria. The exception 
is Sym C which has the same symmetry breaking criteria as 
Sym A, but with twice the coupling strength. Due to this 
we will temporarily ignore Sym C when studying the effects 
of the symmetron parameters. 

We immediately see that the lower assB is, meaning the 
fifth-force has been acting upon the Universe for a longer 
time and also screens at a higher density, the bigger the de¬ 
viations from ACDM are in the extremities of the density 
profile plots, with the differences between the models being 
smallest near the virialization radius. One might naively ex¬ 
pect that models with a lower assB should break off from 
ACDM closer to the centre of the halo, where the density is 
higher, however it is clear from these profiles that the length 
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Figure 2. The top figures show the deviations from ACDM for our f(R) models for CDM and baryons respectively for all density profiles 
with a mass between 10 14 /i _ 1 Mq and 5 X 10 14 k _ 1 Mg. The bottom figures show the same for for all density profiles with a mass 
between 10 13 h _1 Mg and 5 X 10 13 k _ 1 Mg. 


of time that the fifth-force has been working on the particles 
is the stronger of the two effects. 

Returning to Sym C, we note that the increased cou¬ 
pling strength amplifies the amplitude of Sym A, while also 
making it behave as if it had a slightly lower ossb than 
it actually has, in short it simply boosts the effects of the 
fifth-force. Looking at the low-mass figures we note that the 
profiles follow the same trend as before, but that the effects 
from the stronger coupling in Sym C is much stronger. This 
is due to the fact that the halos of this mass range have sizes 
that are approximately of the order of the force range and 
we get a stronger resonance. 

For the f(R) case the main difference between the pa¬ 
rameters is the range of the fifth-force, and the behaviour 
of these f(R )-gravity models have been discussed before in 


Lombriser et al. (2013). What is worth noticing is that we 


observe the same behaviour when comparing the dark mat¬ 
ter density profiles to the gas density profiles as we did in 
the symmetron case. The only exception is for the high mass 
halos in the FofR06 simulation, which is the model closest 
to ACDM, as is expected, in the dark matter case. Flowever, 
the gas profiles show slightly lower clustering in the inner 
regions than what we have for the dark matter (similar to 
what we see for Sym B, but to a larger extent). 


5 TEMPERATURE PROFILES 

In this section we provide the temperature-profiles for the 
halos in our simulations. The motivation for showing this is 
that the temperature of the gas in our Universe is one of 
the direct, and easiest, observables that we can use to com¬ 


pare to our results, furthermore this is something that, to 
our knowledge, has never before been presented for modi¬ 
fied gravity theories. To calculate the temperature from the 
simulations we utilise the ideal gas law, 


p = R a pT , 


(33) 


where p is the thermal pressure, R 3 = -^ s - is the specific 
gas constant and p/p is provided to us from the simulation 
outpuf] In Fig. [ 3 ] and Fig. [ 4 ] we show both the profiles and 
the deviations in the profiles with respect to ACDM for our 
symmetron and /(-R)-gravity models respectively. 

For the models we have simulated the screening prop¬ 
erties for (dark matter) NFW halos have been studied in 
Gronke et al. (20141 and these results are very useful to 


understand the results we find. 

The first thing to notice from the temperature profiles 
is that the effect of modified gravity can be larger than what 
we find in the density profiles and power-spectra. This is a 
similar kind of signature as has been found in the velocity 


field in modified gravity simulations previously (Li et al. 


2013) and not surprising as the temperature of the gas is 


closely related to how fast (and faster means more turbulent) 
the gas is moving. 

The largest halos we study are so massive today that 
they have been screened (as modifications of gravity is a late 
time effect in our models) during most of their evolution. 
This is also reflected in the temperature profiles for FofR05, 
FofR.06 and all the symmetron models where we see a close 
to zero deviation in the centre of the largest halos. The only 


3 


The profiles we show are mass-weighted, i.e. T = 


E ™i T i 

E m i 
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Figure 3. The left figure show the gas temperature profiles for ACDM and our symmetron models for all halos with a mass between 
10 14 h~ 1 MQ and 5 X 10 14 h~ 1 MQ, the right figure show the deviations from standard ACDM for the same models. 
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Figure 4. The left figure show the gas temperature profiles for ACDM and our /( R) models for all halos with a mass between 10 14 
h~ 1 M@ and 5 X 10 14 h ~ 1 Mq , the right figure show the deviations from standard ACDM for the same models. 


exception is FofR04 where we see a deviation even in the 
center. As seen in Fig. 6 of Gronke et al. (20141 all of our 
simulations, except FofR04, have a large degree of screening 
for our largest halo mass range. FofR04 does not have much 
screening even for the largest halos and consequently the 
modifications of gravity are active inside our most massive 
halos leading to a non-zero deviation even in the centre. 


FofR04 is our simulation that is closest to act as a lin¬ 
ear (non-screened) model. For this reason the modification 
in the temperature profiles is close to a constant over the 
whole profile for both mass-ranges. This is to be expected 
since if we neglect the finite interaction range of the fifth- 
force then the modifications of gravity can be thought of as 
just a rescaling of the strength of gravity (i.e. of Newton’s 
constant). For FofR04, the difference in amplitude in the 
temperature with respect to ACDM is seen to be slightly 
larger in the center for our smallest halos. This is likely 
related to the range of the fifth-force, being density depen¬ 
dent and smaller for denser objects, making the effect of 
the fifth-force (through not much screened) smaller for the 
largest halos. 

For FofR05 we see a modification in the centre only for 
our smallest halos. This is also to be expected from Fig. 6 in 


Gronke et al. (2014) which shows that for FofR05 our large 


halos are very much screened, but our small halos are not 
much screened. For FofR06 we see only a very small (a few 
percent) deviation in the temperature for both of our halos 
mass-ranges over the whole profile again consistent with the 
findings of Gronke et al. (2014 1 . The FofR06 model has the 
highest amount of screening and generally the shortest range 
of all models and to see sizeable deviations in this model 
we would need to go to much smaller halo masses which is 
beyond the resolution of our simulations. 


For the symmetron models we also see the effect of the 
additional symmetron screening mechanism. Sym D is the 
symmetron model that has the least amount of (non-linear) 
screening and the fifth-force has been in operation for the 
longest time. However, opposed to FofR04, the deviation in 
the centre of the largest halos is rather small. The reason 
for this is that the coupling to matter )3(<p) oc ip is field- 
dependent and goes to zero for large and dense halos giving 
rise to additional screening which we don’t have the f(R) 
model where /3 is a constant. This effect suppresses the mod¬ 
ifications of gravity inside large halos even though we don’t 
have a large degree of (non-linear) screening. For the other 
symmetron models (Sym A mostly) we also have a large 
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Figure 5. Gas temperature deviations from standard ACDM for 
our /(f?) models for halos with a mass between 1 0 1 3 h~ 1 A/g. and 
5 X 10 13 h~ 1 M Q . 
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Figure 6. Gas temperature deviations from standard ACDM 
for our symmetron models for halos with a mass between 10 13 
1 Mq and 5 X 10 13 k _1 Mg. 


degree of non-linear screening further suppressing the fifth- 
force and thereby suppressing the deviations seen in the tem¬ 
perature profiles. 

The temperature profiles show all the different aspects 
of screening we have in the models we have simulated and 
because of this we would also expect to see an environ¬ 
ment dependence of the temperature profiles for halos re¬ 
siding in different environments similar to what was found 
in Winther, Mota & Li (2012) for dark matter halos. This 


is left for future work. 


6 POWER SPECTRA 


The power spectra were computed using the publicly avail¬ 
able POWMES code (Colombi & Novikov 20111, for both 
the dark matter and the gas. In order to use POWMES to 


compute the gas power spectrum we needed to extract the 



Figure 7. We show the the deviation of the bias from ACDM, 
for f(R) gravity. 



Figure 8. We show the the deviation of the bias from ACDM 
for the symmetron model. 


mass of each cell from the density of the cell. This is done 
by assigning a mass to grid cells using 

m = pVceii. (34) 

where p and Veil correspond to the gas density and volume 
of each cell. 

We will start by studying the effects modified gravity 
have on the ratio between the gas power spectrum and the 
dark matter power spectrum, 

b= (35) 

-t gas 

In Fig. [7] and Fig. [8] we show the fractional difference of 
the gas-dark matter bias with respect to ACDM for both 
the f(R) and symmetron simulations respectively. We first 
note how the bias is strongly scale dependent. At large scales 
(fc < 0.5 h/Mpc) we see that all theories converge to the bias 
of ACDM. However once we move to smaller scales k > 0.5 
h/Mpc we see that the deviations of the bias grow at an 
increasing rate for most models. 

Looking at Sym A, B and D (again, Sym C mirrors Sym 
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k (h/Mpc) k (h/Mpc) 

Figure 9. The left and right figures show the power spectrum deviations from standard ACDM for our various /( fi) models for CDM 
and baryons respectively. 




k (h/Mpc) k (h/Mpc) 

Figure 10. The left and right figures show the power spectrum deviations from standard ACDM for our various symmetron models for 
CDM and baryons respectively. 


A, but is more extreme) we see that the bias can both be 
higher and lower than the ACDM bias depending on ossb- 
The bias seem to rotate in a clockwise motion as ossb in¬ 
creases. We also see a direct link between the strength of 
the coupling and the increased deviations in the bias when 
comparing Sym A and Sym C. For the f(R) simulations we 
see that FofR04 has a bias that increases strongly when we 
go to smaller scales, whereas FofR05 and FofR06 has a bias 
that is close to ACDM. It should be pointed out that Sym A 
and Sym B together with FofR05 and FofR06 are the most 
realistic models, whereas FofR04 and Sym D are extremes 
(which most likely can be ruled out already based on current 
data). The main thing we observe is that the behaviour of 
the bias can be very different depending on how long the 
fifth-force have been in operation, the range over which it 
acts, the coupling strength to matter and the amount of 
screening that is in play. This makes it hard to give a gen¬ 
eral prediction for what the gas does in these theories as it 
is very dependent on the model parameters. However, we do 
notice that it is in the simulations where we have a lower 
degree of screening (Sym D, FofR04 and to a lesser degree 
Sym B) that the relative bias grows significantly with scale. 


A possible explanation for this is that it is only the very 
dense gas inside clusters that experience much screening of 
the fifth-force and therefore do not cluster as expected. 

In Fig. [To] and Fig. i we present the deviations of the 
dark matter and gas power spectra for the symmetron mod¬ 
els and the f(R )~gravity theories from ACDM respectively. 
For both models we see that the majority of the models 
peak at approximately 1 h/ Mpc, which is the scale at which 
streaming velocities reach a value large enough to prevent 
further clustering. The results for dark matter is found to 
be in good agreement with existing simulations in the liter¬ 
ature. 

We start by noting a direct correlation between the am¬ 
plitude of the power spectrum deviations and the time since 
symmetry breaking. The longer time the fifth-force has had 
a chance to work on the dark matter particles and the gas 
the larger deviations from ACDM we find, in agreement with 
expectations. 

Further we see that the higher ossb we have, the lower 
the scale at which the screening mechanism starts, we might 
have expected a complete opposite effect as the density usu¬ 
ally increases as we go to smaller and smaller scale, however 
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this is clearly not what we observe. This might tell us that 
the screening mechanism is less effective for these values 
than naively expected. 

By looking at Sym C we can see the sensitivity of the 
power spectrum with regards to the strength of the fifth- 
force. We observe that the effect of doubling the strength 
of the force compounds drastically in the non-linear regime, 
particularly for the gas component. 

For the f(R) case we know that in the FofR04 simula¬ 
tions we have much less screening than in the FofR05 and 
FofR06 simulations (as can be observed by comparing full 
simulations with linearised no-screening simulations, see|Li,| 


Mota & Barrow (20111). It is only in the most dense re¬ 


gions of the simulation box where the fifth-force is signif¬ 
icantly screened for FofR04. As the most dense regions in 
the simulation box largely consists of highly clustered bary- 
onic gas the gas will experience much more screening than 
the dark matter counterpart and this effect can explain the 
large growth of the bias with scale as we observe for FofR04. 
For models where the screening mechanism works effectively, 
like for FofR05 and FofR06 the fifth-force seem to affect the 
dark matter and gas equally well. 


7 DISCUSSION 

In this paper we have studied the effects modified gravity 
has on structure formation with dark matter and a baryonic 
gas. The aim being to study the effects of modified gravity 
on the baryonic physics in its simplest form. 

We ran several simulations with 256 li particles for the 
symmetron and f(R)~ gravity models respectively. With the 
data from these simulations we analysed the density pro¬ 
files and power spectra for both the dark matter and gas 
components, while also commenting on the behaviour of the 
gas-DM bias. 

We have shown that great care must be taken when 
trying to rule out modified gravity theories by comparing 
observations to the dark matter predictions, as the addi¬ 
tion of the baryonic gas reduces the deviations from ACDM 
significantly, at least for density profiles. The smaller effect 
of modified gravity on the gas density profiles compared to 
the effect it has on the dark matter density profiles is most 
likely due to an environmental effect from the dark matter. 
Namely that the total density p tot = Pdm + Pgas gets suffi¬ 
ciently high to trigger the screening mechanism before the 
gas component has had enough time to be acted upon by 
the extra force to give the expected impact. The collision¬ 
less nature of dark matter give rise to a lower density at the 
inner regions of the dark matter halos than the inner region 
density of the gas halos. 

The analysis of power spectra shows us that it is mainly 
in the non-linear regime that the differences between the 
dark matter and the gas cases are significant, in agreement 
with previous findings. The power spectrum is highly sensi¬ 
tive to the range of the scalar field, and especially sensitive 
to the time that the fifth-force has had to affect the content 
of the Universe. 

We have found that temperature profiles of clusters can 
be a strong signature of modified gravity. The deviation from 
the ACDM predictions can be a factor of a few larger than 
the same deviations found in the density profiles and power 


spectra, and much more sensitive to the fifth-force then pre¬ 
viously assumed. It is also important to notice that the ef¬ 
fects that modified gravity theories have on the baryonic gas 
does not always exactly mirror the effects it has on the dark 
matter. 

Using these newfound characteristics of the symmetron 
and f{R)~ gravity models of modified gravity on the temper¬ 
ature we can start comparing to probes of gravity such as in 


Terukina et al. (20141 or as suggested in Jain et al. (20131 


and hopefully produce better constraints. 

It is worth to note that our study have only focused on a 
simple baryonic gas without taking into account important 
physical effects we know take place in our Universe. In the 
future we aim to extend this work by adding such effects 
as feedback and cooling into our code. Feedback effects will 
particularly be interesting to see whether they enhance the 
deviations we already observe or not, whereas cooling might 
have significant effects on the extremities of the temperature 
profiles, where we observe the extreme deviations. 
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